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Abstract. Weyl's formulation of quantum mechanics opened the possibility of 
studying the dynamics of quantum systems both in infinite-dimensional and finite- 
dimensional systems. Based on Weyl's approach, generalized by Schwinger, a self- 
consistent theoretical framework describing physical systems characterised by a finite- 
dimensional space of states has been created. The used mathematical formalism is 
further developed by adding finite-dimensional versions of some notions and results 
from the continuous case. Discrete versions of the continuous Gaussian functions have 
been defined by using the Jacobi theta functions. We continue the investigation of the 
properties of these finite Gaussians by following the analogy with the continuous case. 
We study the uncertainty relation of finite Gaussian states, the form of the associated 
Wigncr quasi-distribution and the evolution under free-particle and quantum harmonic 
oscillator Hamiltonians. In all cases, a particular emphasis is put on the recovery of 
the known continuous-limit results when the dimension d of the system increases. 
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The continuous Gaussian functions 



9k ■ 



9k{x) = e 2" 



KG (0, 00 ) 



play a fundamental role in physics, particularly in quantum mechanics, due to their 
remarkable properties, among which we mention: 



(i) With the Fourier transform defined as 

poo 



mm 



2,71 



we have 



Jk « 



(ii) The function gi is an eigenfunction of a second-order differential operator 

d 2 



dx 2 



+ x e 2 



e 2 X . 



(1) 



(2) 



(iii) The function g K is a minimum uncertainty state for the coordinate-momentum 

AxAp=^ (3) 

(iv) The Wigner quasi-distribution [25] corresponding to g K defined as 

1 f°° 

W(x, p) = - du e 2ipu g K (x - u) g* K (x + u) 
n J -oc 

is, up to a multiplicative constant, a product of two Gaussian functions 

W(x,p) = 2J^-e- Kx2 e~~y. (4) 
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Figure 1. The functions gi/3 (left), 0x (center) and 33 (right) in the case d = 31. 
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Let d — 2s + 1 G {3, 5, 7, ...} be a fixed odd positive integer, and let Z rf = Z/oK 
be the ring of integers modulo d for which we use {— s, — s + l, s — 1, s} as a set of 
'standard' representatives. The number d represents the dimension of the Hilbert space 
describing the states of the investigated quantum systems. 

The Jacobi theta function [TT1 |2T[ [21] 



has several remarkable properties among which we mention: 

6 3 (z + m + nr, r) = e" iW e" 2 ™ 2 6 3 (z, r) 

3 (z,ir) = -|^exp~^ 9 3 (^-, - 

and [19] 

k in\ 1 v—->. 2ni, fn i 



d 1 d J \f~Kd~ jt^„ \d' nd 

For any k G (0, oo) the function 

g K : Z d — ► R, g K (n) = ^ 



oo 

e 

Q = — OO 

can be expressed in terms of the Jacobi function # 3 as 
and by using the finite Fourier transform 

* n=— s 

Ruzzi's relation ([5]) can be written in a form identical to ([!]), namely, 

F[ Sre ] = ^ 5 i. (6) 



This property of the function $j K and the shape of its graph (see figure [T]) show 
that it can be regarded as a finite version of the Gaussian function g K . We call q k a 
finite Gaussian, and our main purpose is to prove the existence of a finite version for 
the relations (l)-(4). We recover well-known results for continuous Gaussians for large d 
values, i.e. at the discrete-continuous transition, and emphasize conditions under which 
the evolution of continuous Gaussian states can be retrieved from the evolution of their 
discrete counterparts. 

The finite Gaussians g K represent a generalization of Mehta's function f , and the 
relation (|6]) is a generalization of Mehta's relation F[fo] = J'q. More precisely, Mehta 
has proved pj)] that the functions 



00 / 1 — \ 

fk(n) = £ e -%(<*W*H k U% (ad + n)) 

rv= — rv~i > * ' 
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defined by using the Hermite polynomials are eigenvectors of the finite Fourier transform 

F[f k ]=i k f k for any ke{0, 1,2, ...}. 

Since H (x) = 1, we have / = 0i and the relation F[f ] = f coincides to F[qi] = g\. 
Mehta's results are mainly based on the following remarks, which we will use in the 
following: 

a) A periodic function with period d, namely, 

OO f I x 

$:Z-^R, $(n)= J2 vU^fiad + n)) 

a=— oo ^ v ' 

can be defined by starting from any function if : R — > R for which the series is 
absolutely convergent. This is similar to a Zak [26] or Weil [23] transform. 

b) The relation 

°° i>(a+l)\/2nd ^00 

/ </?(x) dx = f(x) dx 

is true for any function if : R — )■ R for which the integral is convergent. 

c) The relation 



is satisfied for any function if : Z — > R for which the series is absolutely convergent. 

In section |2] we review some elements of the mathematical formalism used in the 
case of the quantum systems with finite-dimensional Hilbert space, in a form suitable for 
our purpose. We show that in the case of the free evolution, the time dependent state 
is periodic in time. The behaviour of the commutator of the position and momentum 
operators in the limit of large d is investigated in section |3j Ruzzi has obtained the 
relation ([6]) by using the properties of ^-functions. In section H] we present an elementary 
proof based on a)-c) for this finite version of relation ([!]), and show that Qi is almost 
a minimum uncertainty state. In section |5] we investigate numerically the quantum 
oscillator Hamiltonian. The tendency to have equidistant energy levels becomes evident 
only for d large enough. In section [6] we show, for the first time to our knowledge, that 
the existence of commensurate or equidistant energy levels is a sufficient condition for 
the occurrence of revivals. 

In section [71 by using a)-c) as mathematical tools, we prove that the Wigner function 
corresponding to a finite Gaussian can be written as a sum of four products of finite 
Gaussians. The obtained formula is our main result, and can be regarded as a finite 
version of the relation (jlj). In the particular case k = 1, an expression of the Wigner 
function corresponding to a finite Gaussian Qi has been previously obtained by Marchiolli 
and Ruzzi [12]. They proved that, up to a multiplicative constant, the Wigner function 
corresponding to 
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is 



-tttj/ \ n fn i \ n (2m 2i\ (n i \ n (2m 2i 

where 6*2 and 6*4 are the Jacobi functions 

00 2 
6 2 (z,t) = ^ e i7rr ( a +§) e 27ri ( a +§) 2 

and 

00 

^,r)= £ (-1) 



\Qgi7rTa 2 g27riQ2: 



a=— 00 



2. Quantum systems with finite-dimensional Hilbert space 

In the case of a quantum particle moving along a straight line, the possible positions form 
the set K. = (—00, 00), and the space describing the states of the system is the infinite- 
dimensional Hilbert space L 2 (IR) of all the square integrable functions ip : R — > C. 

We obtain a very simplified version of this continuous one-dimensional system by 
assuming that we can distinguish only a finite number d of positions for our particle. In 
this simplified version, the space describing the states of the system is the ci-dimensional 
Hilbert space C d . Assuming that d is an odd number, d = 2s + l, the space C d can be 
identified with the space % of all the functions 

ip : {— s, — s+1, s — 1, s} — > C 
by using the one-to-one mapping 

H — ► C d : ^ (i>(-s),i>(-s + 1), ...,i[>(s - 
We choose an orthonormal basis {|n)} ne z d in % and define the 'position' operator 

Q-.H^n, Q = sM £ n\n)(n\. (7) 

v n=—s 

The finite Fourier transform 

1 s 

F-.n^H, F = 4= Yl e2fnn '\n){n'\ (8) 

n,n'=—s 

allows us to consider a second orthonormal basis {\k)}kez d , where 

1 s 

\k) = F\k) = F+\ - k) = — V e^ kn \n) 

yd 

v n=—s 

and to define the 'momentum' operator 

P:U—>% p = J?f J2 k \k)(k\. 

k=—s 
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The operators Q and P have the same spectrum, namely, 

/2tt , , /2vr . . /2vr /2tt 



Since 



/27T t / , \ /2vr 

hm W — = and hm (±s) \ — — — ±00 



in the limit d — > 00, the spectra of Q and P correspond in a certain sense to (—00, 00). 
Each state EH can be expanded as 

s s 

\ i))= Y^4>{n)\n)=Y,m\k) 

n=—s k=—s 

where the functions ip:1*d — >C: n tp(n) and ip'.'Ld — >C: k tp(k) satisfying 
tf(n) = (*#> = 4j E ^ kn m , ^ (A) = W> = 4E e-^ fc >(n) . (9) 

are the corresponding 'wavef unctions' in the position and momentum representations 
[5]. The operators Q and P satisfy the relations 

FQF + = P FPF + = -Q. 

The displacement operators A ,B:T-L — > H 



A = e -iV^ p = J2 e-^ k \k)(k\ , B = j^l~ dQ = e^VX^I , (10) 

k=—s l=—s 

are single- valued and satisfy the relations 

A\£) = \£+l) , A\k) =e-^r k \~k) , A d = B d = I , 

B\t) = e^ e \£), B\k) = \k+1), A a B^ =e~^B^A a . 

The general displacements operators j22j EQ] 

D(a, 0) = e^ alB A a B 13 where (a, /3) G Z d x Z d (11) 

define a projective representation of the finite Weyl group. The vectors {|a, 0)} 8 a r=- s i 
where 

\a,fi) = D(a,/3) = V e^« fll (j-a) |j> 



satisfy the resolution of identity[5l 
1 s 

- \<x,P)(a 7 P\ = I. 



d 

a,p=—s 



The tight frame {\a, (3)} s a B=- s can be regarded as a finite system of coherent states [7] 
labeled by using the set Z^xZ^, directly related to the finite phase space Sd><Sd- 

The mathematical objects defined above correspond in the large d limit to those 
usually considered in the case of the quantum harmonic oscillator. An extensive list 
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Table 1. Correspondence between continuous and (2s + l)-dimcnsional case. 



Continuous case 



(2s + l) — dimensional case 



(—00, 00) 



{x\x') = 8{x - x') 
\ip) = J dx ip(x) \x) 

ip{x) = (x\i>) 

x = J dxx \x)(x\ 

x \x) ~ x \x) 
F=^Jfdx'dxe lx ' x \x'){x\ 
= ^=fj dx' dxe lx ' x i){x) \x') 
F\i>) = J dx' T[^}(x')\x') 
F[ip]{x') = -|= / dxe ix ' x ip(x) 
\p)=J r \p) = ^Jdx^ x \x) 
(x\p) = ^e^ 

m')=s{jp-p') 

\i>) = Jdp-ip(p) \p) 
4'(p) = (p|V>) =^jdxe- i ^{x) 
P = J dpp\p)(p\ 
P\f>) =P \P) 



-lap I 



e 2 

= \x+a) 



-lap 



\p) 



-lap 



\i>) 



e i/3x \x) 



e^ x \p) = \p+0) 



$f c dz\z){z\=I 

B2 



V d 

(n\n') = 6 nn > 

M = EIU-.lK»)l») 
ip(n) = (n\ijj) 

Q = "/¥ t n\n){n\ 
Q\n) 



F = J- V s 



nJ 2 -f\n) 



. e a 



|n')H 



^) = ^E^ 1 „=-,e^"'"VW|n') 
F[^](n') = ^E;=- s e^"'"V(n) 



kn I 



m = Cm = ^ d E: = 



■ A;;; 



■0(n) 



p = Jt E k\k)(k\ 

k— — s 



P\k) = kyJ^\k) 
Z-^a— — oo 

A a \n) = \n+a) 
A a \k) = e- 2 Tr ak \k) 



(ad+n) 



= |fc+/3) 
L»(a,/3) = A a B f} 
|a,/3)=£>(a,/3)^|| 
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Example. The Hamiltonian 

H iiee :H^H, H bee = -P 2 
admits the non-degenerated ground level 

Ao = with |0) a corresponding eigenvector 
and the double-degenerated energy levels 



Ai = 


x i2 
d 1 


with 


l±l) 


orthogonal ei{ 


j;envectors 


A 2 = 




with 


|±2) 


orthogonal ei{ 


j;envectors 


A s = 


l* 2 


with 


\±s) 


orthogonal ei{ 


$envectors 



Therefore, 

n=—s 

and the evolution operator 

n=—s 

is periodic with period 2d 

For any state \ip) EH, the corresponding time dependent state 

$:Z d xR — >C, V(n,t) = (n|e _i * Hfrcc |^) 
is periodic in time 

V(n,t + 2d) = y(n,t). 

Note that in the continuum limit, when d tends to infinity, the periodicity of free 
evolution practically disappears, and one obtains the result known from continuous- 
configuration quantum mechanics. 

A similar periodicity has been obtained in [3] for a free wavepacket moving in a 
discrete quantum phase space. However, in [3] the discrete-eigenvalues position and 
momentum operators were defined differently, with the result that the revivals appear 
for minimum uncertainty states but are only approximate for long time evolution in 
other cases. 
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In this section we present a result similar to Floratos [6], in a version adapted to our 
operators Q and P, and a numerical estimation of the spectrum of [Q, P]. The matrices 
of Q and F in the basis {\£)}eez d are 

Since P = FQF + the matrix elements of P in the same basis are 

If we differentiate with respect to x the identity 

Y e fc - = — — 

k=~s 

true for i^O, then we get the relation 

V k e kx = -s e— e ^~ 1 + e— ^(e* - 1) - eV* - 1) 



k=—s 
2m / 



which for x = ^ (j — becomes 



e¥«-0-l" 



2iri 

fc = — S 



Since s = the last relation can be written as 

y ke™w° = (-i) j - £ — — = — 



and we have 

r _o if j 

Pjt-\ _ i /&T (-1)'-' 



2 V d sinf(j-<?) 



if jVi. 



The matrix elements of the commutator [Q,P] are 

f if j=£ 

and for large d, they can be approximated as follows [B] 
[Q,P] je «i(-l)^(fy-l). 

The matrix 

p (-1)^-1));^ 

has the eigenvectors 

/ 1 . an \* i f o r fce{-a,...,-l,l,...,s} 

I -^( — l) J e d J I with eigenvalue < 

W J i=~ s I (l-d)i for fc = 0. 
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Table 2. The eigenvalues r\ n of the commutator [Q,P] in the case d = 15. 



k 


Vk 


k 






k 


Vk 





-27.276466375122 i 


5 


0.999998706977 


i 


10 


1.000016906603 i 


1 


-4.322222514423 i 


6 


0.999999996717 


i 


11 


1.001534631543 i 


2 


0.649632619978 i 


7 


0.999999999998 


i 


12 


1.067898771074 i 


3 


0.988901431861 i 


8 


1.000000000091 


i 


13 


2.560890405316 i 


4 


0.999822475466 i 


9 


1.000000076444 


i 


14 


18.32999286747 i 



This means that d — 1 of the eigenvalues of [Q,-P] are equal to i for large d. We can 
consider that, in a certain sense, 

[Q,P]^i. 

A numerical estimation of the eigenvalues of the commutator [Q, P] in the case d — 15 
can be seen in Table [3 Already for this relative small d value a significant number of 
eigenvalues tend to i. 

4. Minimum uncertainty states 

Let kg (0, oo), sg{1,2, 3, ...} and let d—2s+l. The function 

G K : 1R — > R, G K (x)= e-HV¥M+*)) = £ e -^+*) 2 

Q= — OO Q?=— OO 

obtained by starting from the Gaussian function 
g K :R^R, g K ( x ) = e-f* 2 

is a periodic function with period d. The finite Gaussian 

oo 

g K : Z d — ► E, g„(n) = ^ e"^^ 2 

a=— oo 

which can be regarded as a finite version of g K , satisfies the relation 

0«(-n) = 0«(n), for any n G Z d . 
Theorem 1 [19]. We have 

F[Q K ] = ^gi for any «e(0,oo). ( 13 ) 

Proof. The function G K (x) admits the Fourier expansion 

oo 
l=-oo 

with 
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J- / -2zLife 



' 2tt 
2 V V d 



He/ 



27ri 

e d *""e 



By denoting i = J^f (ad+x) and using the relation 



e s e 



e 4a 



we get [15 



i ™ r(«+i)^ p -^(*V^- ad ) P -ft 2 ^ 



4= e W V d e 2* dt = -He «<* 

/2-n-d J -co Vk5 



whence 



Particularly, we have 



1 oo 

^ E 



e d e Kd . 



e=-oo 



9,(3) =G K (j) = j rd Z7=-oo^ e e 



1 V^s ^-\oo 



27ri 

e d 



-j(ad+n) e --^(ad+n) 



/Kd *-^n=-s A^ct=-oo 



e -^(«H-j) a 



□ 



whence 

.. s oo oo 

* n=—s a=— oo »=— oo 

Since K (— n) = re (n) and Qi(—n) = fli(ra) we have 

s s 

J2n(Q K (n)) 2 = J>(£|i(7i)) 2 = 0. 

n=—s n=—s 

Therefore, the square of the dispersion of Q in the state described by the finite Gaussian 



g K : Z d — ► R, g K (n) = £ 



e _«L M+n) 2 



is 



and, in view of the relation P = FQF + , the square of the dispersion of P is 

27tE: = -^ 2 (0iH) 2 



(ap) 2 = (p 2 ) - (py 



d E:=- S (fli(n)) : 
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Table 3. The state g± is a quasi-minimum uncertainty state. 

d AQAP ||([Q,P])| AQAP-i|([Q,P])| 

3 0.44259776311852 0.44259776311852 0.0 • 10~ 75 



5 


0.49709993841560 


0.49620649757954 


0.000893440 






7 


0.49985914364743 


0.49985140492777 


7.738719663 


io- 


6 


9 


0.49999327972581 


0.49999098992968 


2.289796128 


io- 


6 


11 


0.49999968416091 


0.49999965440967 


2.975123667 


io- 


8 


13 


0.49999998532738 


0.49999998026367 


5.063715121 


io- 


9 


15 


0.49999999932443 


0.49999999924381 


8.061781262 


io- 


11 



We have 



AQAP 




,™ 2 (0«W) 2 



-M*)) 2 \ EU.(Bi(«)) 2 ' 
As concerns the expectation value of [Q, P] , by using the relation (I12p . we obtain 



(Bk\Bk) En=- S (S«H) 2 Vi=-^u=^ sinf(j-f) 
The well-known uncertainty relation originating from Schwarz inequality 

AQAP>±\([Q,P))\ (14) 

is satisfied, and for k — 1 the difference 

AQAP-^\([Q,P})\^0 (15) 

except a few small values of d. Numerical results concerning the case k = 1 are presented 
in table El 

Note that for different definitions of the position and momentum operators in 
a discrete quantum phase space [3], the minimum uncertainty of these operators is 
dependent on the discretization step (the exact formula is known as the generalized 
uncertainty principle) and approaches the result for the continuum case only for an 
infinitely fine discretization; the generalized uncertainty principle can be obtained from 
a quantum mechanical model with discrete eigenvalues for the coordinate operator [2]. 
On the contrary, in our case the uncertainty is approximately minimum for quite small 
d values. Moreover, in [H] the uncertainty relation was shown to reach its minimum 
value only in particular cases, but approximate expansions of the unitary position and 
momentum operators on particular states were used. 



5. Finite-dimensional quantum system of oscillator type 

The Hamiltonian H = ^(P 2 +Q 2 ) is of harmonic oscillator type, but a certain similitude 
between the behaviour of our quantum system with finite-dimensional Hilbert space and 
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Table 4. The eigenvalues of H = \{P 2 + Q 2 ). 



d = 3 


d = 5 


d = 7 


d = 9 


d = 11 


d= 13 


— 


— 


— 


— 


— 


15.685806 


— 


— 


— 


— 


— 


12.088829 


— 


— 


— 


— 


12.908813 


10.202462 










9.802541 


9.713488 








10 156706 


7 964696 


8 21 1687 








7.601849 


7.799516 


7.588461 






7.433857 


5.929737 


6.324626 


6.469345 






5.501405 


5.772956 


5.541025 


5.505452 




4.745031 


4.092770 


4.414645 


4.489404 


4.498956 




3.512928 


3.629951 


3.514121 


3.501381 


3.500114 


2.094395 


2.273277 


2.472337 


2.497725 


2.499837 


2.499989 


1.651797 


1.538153 


1.502561 


1.500166 


1.500009 


1.500000 


0.442597 


0.496978 


0.499856 


0.499993 


0.499999 


0.499999 



the standard harmonic oscillator exists only for a large enough dimension d. Particularly, 
the tendency to have equidistant energy levels becomes evident only for d large enough 
(see Table H] and Figure EJ) • 

The finite Gaussian q 1 is a quasi-eigenstate of the oscillator type Hamiltonian 

H= l -{P 2 + Q 2 ). 



15- ■ 



10 



d=3 d=5 d=7 d=9 

Figure 2. The energy levels of H = \{P 2 + Q 2 ) 



d=ll 



We have (see table EJ) 
Hq 1 w Afli 



for 



A 



HQi (0) 
0i(O) ' 



d=13 



This relation can be regarded as an approximative finite version of the relation 
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Table 5. The state gi is a quasi-eigenstate of H = ^(P 2 + Q 2 )- 







d = 3 


d = 5 


rf = 7 


d = 9 


d = 11 




A 


0.442598 


0.489794 


0.498096 


0.499638 


0.49993 


(H01 


-Afli)(±l) 


2.2-10- 16 


1.2-10~ 2 


2.8-10~ 3 


5.8- 10- 4 


1.1-10" 4 


(H91 


-Afli)(±2) 




-1.2-10~ 2 


-8.5-10" 4 


-1.5-10" 4 


-3.M0~ 5 


(H91 


-Afli)(±3) 






-2.0- 10- 3 


1.3-10" 4 


4.3-10" 5 


(H91 


-A fll )(±4) 








-5.5-10- 4 


-2.9-10~ 5 


(H91 


-Afli)(±5) 










-1.0-10" 4 



6. On the occurrence of revivals 

Theorem 2. If H has k > 2 commensurate energy levels then there exist revivals. 

Proof. If the energy levels £1, £2, Sk are comensurate then £2/^1, £3/21, ••• , £k/ £ i 
are rational numbers and can be represented as fractions. If m is the least common 
multiple of the denominators of these fractions then there exist the integers £2, £3, ... , 
£k such that 

£■ £■ 

-± = -± e Q for any j G {2, 3, k}. 

Ei m 

If l^i), ^2), iV'fe) are the eigenstates corresponding to £1, £ 2 , that is, 

= £j for any j £ {1,2,, fc} 

then the time dependent state 

*:Z d xK — >C, V(n,t) = (n\e- itH \tp) 
corresponding to an arbitrary state of the form (certain coefficients may be 0) 

|V>) = ai\i>i) + "2^2) + ••• + ctk\^k) 
is periodic with the period ^ ZL . Indeed, 

e - itH \tp} = aie" itel |^i) + a 2 e" ite2 |^2) + ••■ + « fc e" itefc |^fc) 



aie" itel + a 2 e- u ™ £l |^2) + ... + a*e~ a ™ £l 



and we have 

e-^+^^He-^l 

whence 

2m7r^ 

Theorem 3. If H has k > 3 equidistant energy levels S\, e 2 , £&, that is , if 



^ ( n,t + ) = 



£2 — £1 — £3 ~ £2 — ••• — £fe — £fc-i 
i/ien t/iere exzsi revivals. 
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Proof. If \?p 2 ), l^fc) are corresponding eigenstates 

Htyj) = Ej \ipj) for any j G {1, 2, , k} 
then the time dependent state 

*:Z d x8 — >C, V(n,t) = (n\e-' ltH \ip) 
corresponding to an arbitrary state of the form (certain coefficients may be 0) 

= ail^i) + "2^2) + ••• + «fc|^fc) 
is periodic with the period ^ ■ Indeed, 

e -' ltH \ip} = aie" itel |^i) + a 2 e- it62 \ip 2 } + ••• + a k e-' lt£k \ip k } 

= e~ it£l («i|^i) + a 2 e-' lt{£2 ~ £l) \ip 2 ) + ... + a k e~' lt( - £k ~ £l) \ip k )) 
= e~ it£l (ai\i) 1 )+a 2 e- it ( £2 - £l ^\ip 2 )+ ... +a k e~ it ^ k - 1 ^~ £ ^\ip k )) 
and, up to a phase factor, we have 

whence 

The number ^ is a period for any coefficients a±, a 2 , ... , a k , but generally, it is not 
a "fundamental" period. For example, in the particular case 

\ip) = a^ifa) + ai 3 \ip 3 ) 

there exists a smaller period, namely, E2 Z ei ■ D 

In the case of our finite-dimensional oscillator, we have equidistant levels and hence 
revivals only for d large enough. Once again, the discrete-continuum transition recovers 
the standard results for the quantum harmonic oscillator. But, more importantly, the 
proposition above relates the revivals with the condition of equidistant energy levels. 
From a physical point of view this relation is extremely important: unlike for free 
evolution, where the revivals period was determined by d but the energy spectrum had 
no equidistant levels and thus the revivals disappeared in the continuous limit, for the 
harmonic oscillator case d does not explicitly enter the expression of the revivals period. 
This period does not disappear in the large d limit; on the contrary, a large d guarantees 
equidistant energy levels, which determine a sort of physical feedback necessary for 
revivals. 

7. Discrete Wigner function 

Let k G (0, 00), sG{l, 2, 3, ...} and let <i = 2s + l. The periodic function : K — > R 

00 2 00 
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with period d allows us to define the function g+ : — > M. 

oo oo 
a=— oo a=— oo 

Since 

0+(n) = G K (n+s+±) and g K (n) = G K (n) 

the function $j+ is a kind of translated finite Gaussian (see Figure EJ). By direct 
computation one can prove the relations 

Qt(~ n ) = at ( n )> 0«( 2ri ) = 04 K H + 0^(n) 

and 



a -if(ad+2n) 2 



04k, (n) - 0^(71) 



Lemma 1. The finite Fourier transform of g£ K satisfies the relation 

1 



n&](2m) 



7 2k 



2(771) - g£(m) 



Proof. The periodic function G^ K (a;) admits the Fourier expansion 

00 



l -tx 



e=-oo 



with 



ct 



f d e-^J2Z^e-<^« a+ ^ d+x )Ydx 



1 v^oo C 

d l^*=-oo Je e d e 



dx. 

By denoting t= \J^j ((a + §) d + x) we get 

= 1 v oo r( a +l/2)V2^d -2fi(t^-(a+l)d 

■ / 2wd ^a=-oo J( a -l/2)\/27rd 



V2nd 

\f2~i~d 
(-l) e roo 



" K * 2 (it 



<ot=— 00 J f 



(a+l/2)v / 27rd 



/2nd J -00 



(a-l/2)V2ird 

d e 



TOO —\(f./£2L — K f 

J___ e V d e Kt (it 



<* e K * (it 



e 2nd 



'2nd 



whence 



-. 00 

= E 



'2nd 



e d 



■lx 



-1)V 



£=-00 



Particularly, we have 

02k CO 



lys y 



E 



00 

£=-co 

12 IT 



e d j*- 



■ l) £ e 2«d^ 2 



00 

a=— 00 



1 V^s 
Vd ^ n = 
1 sr^s 

Vd 2-> n = 



2ni 

e d 
e 



7n (-1)" V- 



2k 



00 

a=— 00 



rE 



2k ^a=-oo 



-l) a e" 



2^ ( ad +") 2 



(16) 
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whence 

(_-\\n 00 

nji.](") = L #- E (-i)°e-^^) 2 

a=— 00 

and we get 

JT^J(2m) = -= E (-l)"e-^^ +2m ) 2 = -= aH-fltH • □ 

Lemma 2. // the numbers N a ^ are such that the series are absolutely convergent then 

00 00 00 

E Na 'P = E N »+V,V-V + E N V+r)+hV-rr 
a,(3=—oo ^,rj=~oo fi,r]=—oo 

Proof. We separate the sum as 



E N ^= E ^+ E ^ 

a,/3=-oo a,/3 

both even one even 

or and 

both odd other odd 

and use the substitutions (a, j3) = (fi + r], fi — rj) and (a,/3) = (// + 77 + 1, \i — rj), 
respectively. □ 

The function W : Z d x Z d — ► C 

W(n,m)=i E e^ mk g K (n - k) g K (n + k) 

k= s~ S 4 . , , (17) 

— I V e^T mfc E e -ftf (ad+n-fc) 2 e -Kf (/3d+n+fc) 2 
k=—s a,/B=— 00 

is called the discrete Wigner function corresponding to K . It is well-determined by its 
restriction to the unit cell {— s, —s+1, s— 1, s}x{— s, — s+1, s— 1, s} directly related 
to the finite phase space SjX-Sd- 

Theorem 4. TTie discrete Wigner function W is a sum of products of finite Gaussians 



+73^0^l r 

Proof. By using theorem 1, lemma 1 and lemma 2 we get 

s 00 

W(n,m) — - J2 e^r mk E e -«f (Gu+^d+n-fc) 2 e - K |(( M -r,)d+n+fc) 2 

fc=— s /i,r]=—oo 
s 00 

+ 1 e ^ imfc E e -Kf ((At+77+l)d+n-fc) 2 e - K f ((/j-^d+n+fc) 2 

k=~s ii,rj=—ca 

* . . 00 
4tti 



18) 



2«f (/i.d+n) 2 p -2 K f { V d~ky 

k=— s /i,ri=—oo 



\ E ^ mk E e «nv— v e ~* 



d 

k=—s n,r)=—ca 
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02k (™) 



02k W 



12 (m)+g1j (m) 



L -.1 



r 



12 (m)-g^ (m) 



Figure 3. The functions involved in relation (p~8]> in the case k = I, d=31 



7S E e 



/i=—oo 
oo 



£ e 

?7=— oo 



-2K^(r)d-k) 2 



fi=—oo 



sfd 



r]=—oo 



k=—s 



r)=—oo 



k=d 2 -( n ) £»^ 2m ) + 73 ^[S^]( 2m ) 



/2kcZ 



'2nd 



□ 



The shape of the functions involved in the expression (1TB]) of can be seen in Figure [3j 
In the finite phase space (which is a finite torus), the obtained discrete Wigner function 
W has three peaks placed around (0,0), (0,s), (s, 0) and an anti-peak around (s, s) 
(see Figure H]). Note the similarity of the representation in Figure H] and the results in 

pjnsiiizi. 

The Wigner function from the continuous case is, in a certain sense, the limit of the 
discrete Wigner function. Therefore, in the continuous limit, the only "surviving" peak is 
expected to be that around (0, 0), in agreement with the known results from the Wigner 
function of continuous Gaussians. The disappearance of the peaks placed around (0, s), 
(s, 0) and of the anti-peak around (s, s) in the continuous limit is rather mysterious. The 
numerical simulations show that for large d all of them keep an amplitude comparable to 
that of the peak placed around (0,0). A possible explanation is the following. In finite 
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Figure 4. The Wigner function W(n,m) in the case k= |, d=31. 

dimension these peaks are placed on the "boundary" of the finite phase space Sd x Sd- 
In the continuous limit they must be located on the "boundary" of the phase plane M 2 
which contains just one point, namely oo. It is known that the compactification of the 
real or complex plane is obtained by adding one point, and the extended plane (called 
Riemann sphere) corresponds to the unit sphere through the stereographic projection. 
We can consider that, for d — > oo, that is, for s — > oo, the peaks placed around (0, s), 
(s, 0) compensate with the anti-peak around (s, s). 

8. Discussions and Conclusions 

The simpler way to define a Gaussian type function on the set {— s, — s + 1, ... , s — 1, s} 
or Sd is to consider the restriction of the Gaussian g K (x) =e~^ x to these sets 

/ : {-s, -s + 1, ... , s-l, s} — ► R, f(n) = e"f" 2 

respectively, 

f:Sd^R, f(n^pJ=e- K > 2 . 

The functions obtained in this way behave similar to the continuous Gaussians only for 
large values of d = 2s + l. 

Our approach is based on an alternative method. By using a Weil-Zak type 
transform, we firstly generate a periodic function G K with period d, and then we define 
our finite Gaussian as a restriction of G K to Z, namely g(n) = G K (n). Our finite Gaussians 
starts to behave similar to a continuous Gaussian from relatively small values of d. In 
addition, they have several remarkable mathematical properties. 

The discreteness of the configuration space has a particular appeal in both classical 
and quantum physics, particularly in phase space p2] or in attempts to unify general 
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relativity and quantum mechanics [12]. On the other hand, the results in a discrete 
configuration space are required to match standard physical results in the continuum 
limit. Because discrete physical systems are of interest in quantum mechanics or in 
classical physics in connection to the problem of sampling, the characterization of such 
systems and of their evolution has already received some attention. For example, a 
discrete quantum phase space has been studied in (81 |£j, the continuous limit being 
obtained by decreasing the lattice spacing in both momentum and position space, while 
in [2] the dynamics on a discrete quantum phase space has been investigated by defining a 
Hamiltonian with the appropriate classical limit. The results in this paper are different 
from previous ones and demonstrate the importance of the choice of some particular 
states in the discrete-continuum transition. As an additional example in this respect, 
we mention that a discrete model of the quantum harmonic oscillator, for instance, 
can be introduced such that the energy spectrum is equally-spaced and the spectra of 
both momentum and position operators are denumerable non-degenerate pQ. This last 
model, too, recovers the results for the ordinary harmonic oscillator in an appropriate 
limit. As do the discrete models of the quantum harmonic oscillator in terms of Kravchuk 
polynomials pU] or Harper functions [3]. 

Summarizing, the correct continuous limit of discrete models can be obtained in 
many situations. The choice of Gaussian states render the discrete- continuous transition 
smoother in the sense that many results known from the continuous case are obtained 
for reasonably large d values or can easily be extrapolated from the results for finite 
d values. Moreover, Gaussian states have advantages over other states in phase space 
representations of physical systems since, as shown in Section 7, the Wigner distribution 
function has a particularly simple expression in this case. In conclusion, finite Gaussian 
states represent a useful mathematical tool for the study of quantum or classical physical 
systems. 
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